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Abstract 

The symmetries, especially those related to the R-transformation, of the reflection equation(RE) 
for two-component systems are analyzed. The classification of solutions to the RE for eight-, six- 
and seven-vertex type R-matrices is given. All solutions can be obtained from those corresponding to 
the standard R-matrices by A'-transformation. For the free-Fermion models, the boundary matrices 
have property tr R+(0) = 0, and the free-Fermion type R-matrix with the same symmetry as that 
of Baxter type corresponds to the same form of AT-matrix for the Baxter type. We present the 
Hamiltonians for the open spin systems connected with our solutions. In particular, the boundary 
Hamiltonian of seven-vertex models was obtained with a generalization to the Sklyanin’s formalism. 


1 Introduction 

In the framework of quantum inverse scattering method(QISM)[Q, ||, ||, [| |[, Yang-Baxter equation(YBE) 

Rl2(u)R 13 {u + v)R -2 3 (v) = R-23{v)Rl3(u + v)R 12 {u) (1.1) 

is a sufficient condition for the integrability of systems with periodic boundary condition(BC). Given 
a solution R-matrix to YBE Q), we can construct the Lax operator of certain models at suitable 
representation of R, and hence transfer matrix t{u). The YBE ensures that t(u ) commutes with each 
other for different spectrum parameters. So, if we expand t{u) with respect to spectrum parameter u, 
the coefficients are a set of conserved quantities which satisfy Liouville’s criterion of integrability § 0 - 

However, when considering systems on a finite interval with independent boundary conditions at each 
end, we have to introduce reflection matrices K±(u) to describe such boundary conditions. Sklyanin 
assumed that R-matrix has the following symmetries^], 

Regularity: R(0) oc P; 

P-symmetry: Pi 2 Ri 2 (w)Pi 2 = R 2 i(w) = Ri 2 (w); 

T-symmetry: R^* 2 = Ri2(w)> 

Unitarity: Ri 2 (u)R 2 i(— u) oc id; 

Crossing Unitarity: R^^R^—it — 2p) oc id, 
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where 77 is crossing parameter and P 12 is the permutation matrix. t\ , £2 denote transposition in space 
Vi and V 2 respectively. In order that the BC are compatible with integrability, the reflection matrices 
should obey so-called reflection equations(RE), or boundary Yang-Baxter equations(BYBE)|| ||, 0 

1 2 

Ri 2 (u - v) K- (u)R 2 i(u + v ) K- (v) = 

K- (v)Ri 2 (u + v) K- (u)R 2 i(u - v), (1.2) 

1 2 

R\ 2 (—u + v) K+ {u)R\ 2 (— u — v — 2 77 ) K+ (v) = 

2 1 

K+ (v)Ri 2 {—u — v — 2rj) I\+ (u)Ri 2 (—u + v). (1.3) 

where K±= K± Cg) 1, K±= 1 ® K± , and R(u) satisfies YBE 0). So for a solution A'_ (u) to RE ( |L2| ) , 
the relation 

K + (u) = K^_(—u — 77 ) (1.4) 

gives the solution to equation ©• Nevertheless, not all P-matrices possesses the above-mentioned 
properties, some generalizations should be made (see e.g.|jlT)). As we will see in section 4, the seven- 
vertex(7V) models are also beyond Sklyanin’s formalism for their P-matrices do not enjoy T-symmetry. 
It was stated in ||l2| that if P-matrix has regularity, unitarity, and crossing unitarity symmetries, but 
does not have PT-invariance, we can propose K + (u) to satisfy the equation: 

1 2 

Pi 2 (—u + v) K + (u)R 2 i(—u — v — 2rj) K + ( v ) = 

I< + (v)Ri 2 (—u — v - 2rf) K + (u)R 12 {—u + v), (1.5) 

and the integrability can be proved as well. There is also an relation between the solutions of Eq. (|l.2| ) 
and Eq.(p~5|) 

K + (u) = K-(—u — 77 ). (1.6) 

We will find later that the Baxter type and free-fermion type I solutions of seven-vertex models are in 
this case. 

Due to the significance of the reflection equations, a lot of work have been directed to the study 
of their solutions 00 Q 0 . BEE- and the Hamiltonians of the systems with such boundary 
conditions are also constructed. However, most of those work are based on the P-matrices which are 
derived directly from the parametrization of the statistical weight in vertex models. There in fact exist 
many kinds of P-matrices according to the classification of eight- and six-vertex type solutions of both the 
YBE and the coloured YBE [|If| ^o[ [|i[] . It is tedious to solve the reflection equation for every P-matrix. 
Fortunately, all those P-matrices of two-component systems can be obtained by applying particular 
solution transformation to standard(or called gauge) ones J 20 ) which satisfy certain initial conditions (let 
us call solution transformation of P-matrix as R-transformation for brevity). The word ’’two-component” 
means that there exist two states in the systems: particles and antiparticles in Held theory, spin up and 
down in spin system, arrow up and down or right and left in lattice model (see Ref. |22|). After a detailed 
study of reflection equation, we can show that there exists a corresponding transformation to the K- 
matrix (we call it A^-transformation) to keep RE invariant under P-transformation. Therefore, we only 
need to concentrate on Af-matrices for the standard P-matrices. 

In this paper we shall focus our attention on solutions to reflection equations of two-component 
systems up to AT-transformation. The solutions are divided into three cases, each of corresponds to 
eight-vertex( 8 V), six-vertex( 6 V) and seven-vertex(7V) models and will be discussed in section 2, 3, 4, 
respectively. For each case, we first analyze the symmetries of the RE, especially those related to the 
P-transformations, and then find solutions to the RE for the Baxter type and free-Fermion type standard 
P-matrices, respectively. We put emphasis on new solutions, but for completeness, we also give the 
solutions obtained by others. 

In section 5, for those solutions given in previous sections, we shall construct the corresponding 
local Hamiltonian of the open spin-chain. The local Hamiltonian means that it only consists of nearest- 
neighbor interaction terms. A system with such Hamiltonian can be viewed as having coupling with 
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magnetic field on its ends. Finally we shall argue that for all boundary conditions to the free-fermion 
models, the reflection matrices K + (u) have property trK + ( 0 ) = 0 . This property requires us to derive the 
Hamiltonian from the second derivative of the transfer matrix |h|. In section 6, we make some remarks 
and discussions. 


2 /\_-Matrices to Eight-Vertex Model 


In this section, we shall first study the symmetries of RE and give the A'-transformation corresponding 
to ^-transformation in [Q. With these discussions, we can concentrate our attention on the RE for 
the standard 8V A-matrices, which are divided into three types: Baxter-type (or XYZ spin-chain [|J), 
Free-Fermion type I (or XY model| 23 |, ^ 5 |), and Free-Fermion type II. All AT-matrices associated with 
these A-matrices are given. 


2.1 Symmetries of Reflection Equation 


The general eight-vertex A-matrix and the corresponding K_(u) matrix are expressed in the following 
forms respectively, 

( 0 0 UJ7(u) ^ 

0 U>2 (u) Ul 8 {u) 0 

0 u>s(u) u>3 (u) 0 

\ U)8(u) 0 0 u>4(u) ) 


R(u) = 


( 2 . 1 ) 


a2 <"{V 

V a 3 \ u ) a 4 :{u) J 


( 2 . 2 ) 


Assuming that A(u) is a solution to YBE (0, then as studied in |20| , there are four symmetries for 
eight-vertex type A-matrix ©■ 


(R.A) Symmetry of interchanging indices. If we exchange the elements of R(u) as aq(it) <-> u>4{u), 
W5 (u) <-> loq(u) or u)2(u) <-> (u), uJr(u) <-> the new matrix also satisfies YBE Q)- 

(R.B) The scaling symmetry. Multiplication of R(u) by an arbitrary function f(u) is still a solution 
to YBE (O). 


(R.C) Symmetry of spectral parameter. If we take a new spectral parameter u = Xu , where A is 
an constant complex number, the new matrix A(w) is still a solution to YBE ( 1 . 1 ). 


(R.D) Symmetry of weight functions. If we replace weight functions u>t(u), ojs(u) by the new ones 

U> 7 (u) = S _1 ui 7 (u), Lu s {u) = sw 8 (u), (2.3) 


where s is a non-zero complex constant, the new matrix is still a solution to YBE ( 1 . 1 ). 


The symmetri es ( R.A)-(R.D) are called solution transformations (or A-transformation) of 8V type 
solution of YBE (|l.l|). 


It is convenient for later discussion to use such notation as follows 


LOi(u — v)= Ui, 0Ji(u + V) = Vi, 
a,i(u) = Xi, Ui(v) = yi. 


Substituting the matrices A and AY into the reflection equation (E 2 ), we get sixteen component equa¬ 
tions, which are divided into groups according to symmetries of the indices: 
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(u 7 v 8 - u 8 v 7 )x 4 y 4 - u 8 v 2 x 2 y 2 + u 7 V3X 3 y3 + u 4 v 6 (x 2 y 3 - x 3 y 2 ) = 0 , 

(■ u 7 v 8 - u 8 v 7 )x 4 yi - u 8 v 2 x 2 y 2 + u 7 v 3 x 3 y3 + Uiv 8 (x 2 y3 - x 3 y 2 ) = 0 , 

(u 2 v 3 - U3V 2 )x 1 y 4 + u 2 v 8 x 2 y 2 - u 3 v 7 x 3 y3 + u 5 v 1 (x 3 y 2 - x 2 y 3 ) = 0 , 

(u 2 v 3 - U3V 2 )x 4 y 1 + u 2 v 8 x 2 y 2 - u 3 v 7 x 3 y3 + u 6 v 4 (x 3 y 2 - x 2 y 3 ) = 0 , 

u 7 V!Xiyi - u 7 v 4 x 4 y 4 - u 4 v 7 x 4 y 4 + u 4 v 7 x 4 y 4 + (u 4 - u 4 )v 2 x 2 y 2 + u 7 (v 5 - v 6 )x 3 y 2 = 0 , 
u 8 vixiyi - u 8 v 4 x 4 y 4 - u 4 v 8 x 4 y 4 + u 4 v 8 x 4 y 4 + (u 4 - Ui)v 3 x 3 y3 + u 8 (v 5 - v 6 )x 2 y 3 = 0 , 

U 3 VeXiyi - u 3 V5X 4 y 4 - u 8 V3X 4 y 4 + usV3X 4 yi + (its - u 8 )v 8 x 2 y 2 + U3(v 4 - v 4 )x 3 y 2 = 0, 
u 2 vexiyi - u 2 v 8 x 4 y 4 - u 8 v 2 x 4 y 4 + u 8 v 2 x 4 y\ + (its - vs)v 7 X3y3 + u 2 (v 4 - v\)x 2 y 3 = 0, 

(itii>i - u 3 v 2 )x 4 y 2 + (u 7 v 8 - u 5 v 5 )x 4 y 2 - u 5 v 4 x 2 yi + u 4 v 5 x 2 y 4 - u 3 v 7 X3y 4 + u 7 v 3 x 3 y 4 = 0 , 

- u 2 v 3 )xry3 + (u 8 v 7 - u 5 v 5 )x 4 y 3 - u 5 vix 3 y 4 + u 4 v 5 x 3 y 4 - u 2 v 8 x 2 y 4 + u 8 v 2 x 2 y 4 = 0 , 

(u 4 v 4 - u 3 v 2 )x 4 y 2 + ( u 7 v 8 - u 6 Vo)xiy 2 - u e v 4 x 2 y 4 + u 4 v 6 x 2 y 4 - u 3 v 7 x 3 y 4 + u 7 v 3 X3y 4 = 0 , 

(u 4 v 4 - u 2 v 3 )x 4 y3 + (u 8 v 7 - u 6 v 6 )x 1 y 3 - u e v 4 x 3 y 4 + u 4 v 6 x 3 y 1 - u 2 v 8 x 2 y 4 + u 8 v 2 x 2 y 4 = 0 , 

uev 2 x 4 y 2 - u\v 7 xiy 3 + u 2 v 5 x 4 y 2 - u 7 v 4 x 4 y 3 + {u 2 v 4 - u 4 v 2 )x 2 yi + (u 6 v 7 - u 7 v 6 )x 3 yi = 0 , 
uev 3 x 4 y3 - u 4 v 8 x 4 y 2 + u 2 v 5 x 4 y 3 - u 8 v 4 x 4 y 2 + (u 3 v 4 - u 4 v 3 )x3y 4 + (u 6 v 8 - u 8 v 6 )x 2 y 4 = 0 , 
U5V 2 x 4 y 2 - u 4 v 7 x 4 y 3 + u 2 v 8 x 4 y 2 - u 7 v 4 xiy3 + (u 2 v 4 - u 4 v 2 )x 2 y 4 + {u 8 v 7 - u 7 v 8 )x3y 4 = 0 , 
U5V3X 4 y3 - u 4 v 8 x 4 y 2 + 113^6^12/3 - u 8 v 4 xiy 2 + (113^4 - u 4 V3)x 8 y 4 + {u 8 v 8 - u 8 v 8 )x 2 y 4 = 0 . 


(Al) 


(A.2) 


(A. 3) 


(A 4 ) 


(AS) 


(A6) 


After a careful study of the above equations, we find that if one apply the following transformations 
to K_(u) under the transformations (R.A)-(R.D), the system of equations (A) keeps invariant: 


(K.A) The symmetry of interchanging indices. This symmetry will be discussed for each type of 
R-matrix later. 

(K.B) The scalar symmetry. If we multiply K_(u) by an arbitrary function g(u), the new matrix 
g(u)K-(u) is still a solution to RE. On the other hand, all the R-matrices up to an arbitrary scalar 
function have the same reflection matrix. 

(K.C) The symmetry of spectral parameters. If we take a new spectral parameter u = Xu , where 
A is any constant, the new matrix K(u) also satisfies RE for R(u). 

(K.D) The symmetry of weight function. If applying the transformation (R.D) to R(it), we can 
make a corresponding A-transformation on R'_(it): 

a 3 (it) = \Zla 3 (u), a 2 (u) = yfs 1 a 2 ( , u) (2.4) 

keeping ai(u), a 4 (u) unchanged. The new K_(u ) matrix is a also solution to RE for new R-matrix. 


Considering the above symmetries for R-matrix and AA-matrix, we can concentrate on the standard 
R-matrix with the restrictions [po| 

w 5 (it) = w 6 (it) = 1, lu 7 (u) = w 8 (w), (2.5) 

and initial condition 

Ri 2 (0)=Pi2. (2.6) 


Note that from condition (2J5), we only need consider R-transformation (R.A) of interchanging indices 
to 1 <-> x> 4 and lo 2 <-> to 3 hereafter. All R-matrices are classified into two classes: Baxter type and free- 
fermion type, according to whether or not the elements of R-matrix satisfy the free-fermion condition 

HH.HI 

u! 4 (u)uj 4 (u) + u> 2 (u)u> 3 (u) — u> 8 (u)u] 8 {u) — u> 7 {u)lo 8 (u) = 0. (2-7) 


The RE corresponding to these two kinds of R-matrices has very different properties. We shall discuss 
solutions to RE for these gauge R-matrices respectively. 
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2.2 Baxter Type 


The gauge i?-matrix of Baxter-type was first derived by Baxter it has the following paramatrization, 

{ oii(o) = 014(0) = sn(zi + h)/ sn h, 

uj 2 {u) = v 3 (u) = sn u/ sn h, , . 

w 5 («) = w 6 («) = 1, 1 } 

u> 7 (o) = oi§(o) = k sn u sn(o + h), 


where snzt, cnzi, dnu are Jacobi elliptic functions of modulus k. It is a high-symmetric one with oii(zz) = 
014(zi), w 2 {u) = w 3 (u), so the transformation of interchanging indices (R.A) has no effect in this case. 

The iv_-matrices in this case have been widely discussed in 610 |l 5 [ |l8|. For completeness, we list 
here main results. The most general one was given in @0 as follows 


K _ ( u ) 


/i sn(2o) 


v sn(a — u) 

A(1 — k sn 2 u) —1 — k sn 2 u 
I - I / 2 sn^a sn’^n 


y sn(2o) 


A(1 — k sn 2 u)-f-l-l-fc sn 2 ii 
1 — k 2 sn 2 a sn 2 v/ 


v sn(a + u) 


( 2 . 9 ) 


where /i, v , A, a are free parameters, and the other special solutions can be obtained by setting these 
parameters to take special values. 


2.3 Free-Fermion Type I 

The R-matrix of free fermion type-I is less symmetric than that of Baxter type. In this case, w 2 (u) = 
w 3 (u), but w±(u ) 7^ u>4 (o). The reflection equation is equivalent to five systems of equations: 


! u 2 v 7 (x 2 y 2 - x 3 y 3 ) + u 5 Vi(x 3 y 2 - x 2 y 3 ) = 0 , 
u 2 v 7 [x 2 y 2 - x 3 y 3 ) + 05114(2:31/2 - x 2 y 3 ) = 0 , 
u 7 v 2 (x 2 y 2 - x 3 y 3 ) + 04115(2:31/2 - x 2 y 3 ) = 0 , 
u 7 v 2 (x 2 y 2 - x 3 y 3 ) + 04115(2:31/2 - x 2 y 3 ) = 0, 

f u 7 v 1 x 1 y 1 - u 1 v 7 x 1 y 4 + U 4 V 7 X 4 yi - 07042:41/4 + (o 4 - ui)v 2 x 2 y 2 = 0 , 
\ u 7 vxxiyi - u\v 7 xiy 4 + U 4 V 7 X 4 yi - 07042:42/4 + («4 - ui)v 2 x 3 y 3 = 0 , 


(B.l) 


(B. 2 ) 


f u 2 v 3 {xxyi - X 4 y 4 ) + 0502(2:41/1 - x\ j/ 4 ) + 02(04 - Vi)x 2 y 3 = 0 , 
\ u 2 v 3 [x\y\ - 2:41/4) + 0502(2:41/1 - 2:12/4) + 02(04 - v\)x 3 y 2 = 0 , 


f u 5 v 2 xiy 2 - u\v 7 x\y 3 + 02052:41/2 
u 3 v 2 x\y 3 - u\v 7 x\y 2 + 02052:41/3 
u 5 v 2 X4y2 - U4V 7 X4y 3 + u 2 v 5 xxy 2 
u 5 v 2 X4y 3 - U4V 7 X4y 2 + u 2 v 5 xiy 3 


u 7 V4X4y 3 + (u 2 vi - uiv 2 )x 2 yi + (0507 - 0705)2-31/1 = 0 , 
u 7 V4X4y2 + ( u 2 vi - uiv 2 )x 3 y\ + (0507 - 0705)221/1 = 0 , 
u 7 v\X\y 3 + {u 2 v4 - u 4 v 2 )x 2 y 4 + (o 5 o 7 - u 7 v 5 )x 3 y4 = 0, 
u 7 viXiy 2 + (o 2 o 4 - U4V 2 )x 3 y4 + (o 5 o7 - 0705)221/4 = 0, 


(B. 4 ) 


! (uiVi - u 2 v 2 )xiy 2 + (u 7 v 7 - U 5 v 5 )x 4 y 2 - u 5 vxx 2 yi + O1O522I/4 
[uivx - u 2 v 2 )x 1 y 3 + (0707 - u 5 v 5 )x 4 y 3 - u 5 vix 3 yx + Uxv 5 x 3 y 4 
(u 4 v 4 - u 2 v 2 )x 4 y 2 + (u 7 v 7 - u 5 v 5 )xxy 2 - u 5 V4X 2 y 4 + U4V 5 x 2 yx 
(u 4 v 4 - u 2 v 2 )x 4 y 3 + (u 7 v 7 - u 5 v 5 )xxy 3 - u 5 v 4 x 3 y 4 + U4V 5 x 3 yx 


U2V 7 x 3 yi + u 7 v 2 x 3 y 4 = 0 , 
u 2 v 7 x 2 yi + u 7 v 2 x 2 y 4 = 0 , , , 

u 2 v 7 x 3 y4 + u 7 v 2 x 3 yi = 0, ^ 
u 2 v 7 x 2 y 4 + u 7 v 2 x 2 yi = 0. 


There also exist symmetries of interchanging indices. The system of equations (B) is invariant under 
exchange of ai(u) <-> 04(11) and u>i(u) <-» 014(11) or a 2 {u) a 3 {u). The gauge I?-matrix is pffl, 

u)i(u) = cn u + H snzt dnzi, 

014(11) = cnu — H snzt dnzi, 

012(11) = 013(11) = G snzt dnzi, ( 2 - 10 ) 

015(11) = 015(11) = dnzi, 

017(0) = oig(o) = k snzt cno, 
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where G, H are arbitrary parameters with relation G 2 — H 2 = 1 . Note that in ( 2 . 10 ) we do not take 
u>5 (u) = loq (12) = 1 in oder to compare our following discussion with other’s work. We will consider the 
general i?-matrix which has wi(u) 7^ 014(12), i.e. H 7^ 0 . The case of H = 0 is remarked at the end of this 
subsection. Now we solve the RE (B) case by case. 

Case 2 . 3 . 1 : Diagonal solution. From (B.l), 02 (u) EOtt <23(14) = 0 . There are only two equations to 
be considered, 

u 2 v 5 (xiyi - £42/4) + u 5 v 2 {x 4 y 1 - X12/4) = 0 , ^ 

UiViXiyi — UiVyXiy4 + I24II7X4I/1 — 1271)4X41/4 = 0 . ' 

Introducing new variable g{u) = ai(14)704(14) and solving g(u) from the above equations, we get a solution 


K_(u) = 


cn u dnu ± ik! snw 0 

0 cm2 dnu =F ik' snii 


( 2 . 12 ) 


where k' is the complementary modulus of elliptic function. Note that the diagonal solution of 8V 
free-fermion type I has no free parameter, which is different from that of Baxter type. 

Case 2 . 3 . 2 : Skew-diagonal solution. If 02(12) ^ 0 , we conclude from (B.l) that 

a 2 (u) = ea 3 (u), e = ± 1 . ( 2 - 13 ) 


Taking 01(14) = 0 , we get from (B. 4 ) 


£4(1221)51/2 - 1271141/3) = 0 . 


With the help of ( { 2 . 13 ) and ( 2 . 1 C| ). the above equation calls for 04(12) = 0 . However, this is contradictory 
to 02(14) ^ 0 as seen from (B. 2 ). So the skew-diagonal solution does not exist due to less symmetry of 
R-inatrix. 


Gose 2 . 3 . 3 : Genera l solu tion. Because (121 — 124)112 = (ui — 114)122, the following equation is obtained 
from (B. 2 ), (B. 3 ) and (|T 3 |) 


127111X11/1 - 121117X11/4 + 144117X41/1 - 127D4X41/4 
-e{u2V 5 (xxyi - X41/4) + u 5 v 2 (x 4 yi - X11/4)} = 0 . ( 2 - 14 ) 


Differentiating the above equation with respect to v and setting 11 = 0 , we can express ai(12), 04(12) as 


where 


J 01 (12) = (F(u) cnu dnu — G(u) snu)p(u) / 2 , 
( 04(12) = ( F(u ) cnu dnu + G(u ) sni4)p(i4)/2, 


fc((l - ekG)d + Hc 2 ) 2 
F(u) = ci H- ^ -sn 12, 

fc(( 1 - ekG)c 2 - fc' 2 ffci) 2 

E{u) = c 2 + —- -pr=~. -— sn 2 12. 

eG — k 


( 2 . 15 ) 

( 2 . 16 ) 


and p(u) is a meromorphic function to be determined. Substituting the above expressions into (B. 5 ), we 
get 

02(12) 

——— = /4 snu cnu dm2, 
p{u) 

and an additional restriction between ci and C2 from (B. 2 ) 

fc ,2 c 2 + Cj = 2 /2 2 (G — ek)/k. 


Therefore, the most general solution is 


I\- (12) 


F(u) cn(12) dn(it) + E{u) sn(it) 
2e /i sn(u) cn(n) dn(u) 


2/4 sn(it) cn(u) dn(ii) 

F(u) cn(12) dn(u) — E(u) sn(12) 


( 2 . 17 ) 
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The result in 0 is a specific case of /i 
be obtained by setting fj, = 0. 


= 1 . It is also easy to find that the diagonal solution ( 2 . 12 ) can 


Remark 2.3.1. There are in fact various parameterizaions of free fermionic 8V / 7 -nratrix, one of 
which is given in | 24 | 

uji(u) = 1 — e(u)e(hi)e(h 2 ), 
u>4(u) = e(u) - e(hi)e(h 2 ), 
uj 2 (u) = e{h 2 ) - e(u)e(hi), 

u 3 {u) = e{hi) - e(u)e{h 2 ), ( ■ 

u 5 (u) = ojq{u) = ^e(hi) sn(hi)e(h 2 ) sn(h 2 )(l - e{u))/ sn(^) 
u 7 (u) = u 3 (u) = —i k^/e(hi) sn(hi)e(h 2 ) sn(/i2)(l + e(-u)) sn(lj), 

where h\ and h 2 are colour parameters, and e(u) is the elliptic exponential: 


e(u ) = cn(it) + i sn(it). 


If we make transformation (R.B) to ( 2 . 1 S| ) with factor function 

\/e(hi)e(h 2 ) sn/ii snh 2 -—^777^ 
snu/2 


and set 


hi = h 2 = h, u —> u/ 2 , G = 


H = 


cnh 


sn h ’ sn h ’ 

the new / 7 -matrix coincides with ( 2.10| ), so our solution includes the diagonal solution given in JTc| . 
Remark 2.3.2. Let us consider the special case of H = 0 and G = 1 . In this case, the / 7 -matrix is 


! u>i(u) = 1x4 (n) = cnu, 
lu 2 (u) = cu 3 (u) = snu dnu, 
W5 (u) = u>6 (u) = dnw, 
to 7 {u) = cos(u) = k snu cn u. 


( 2 . 19 ) 


It has the same symmetry as Baxter type. The calculation shows that this feature is responsible for the 
fact that both two / 7 -matrices share the same K-(u ) as given in (p. 9 [). 


2.4 Free-Fermion Type II 

This kind of / 7 -matrix takes the form, 




lui(u) 

u 2 {u) 


OJ4 (u) = 
-W 3 (w) 


cosh(AM) 
cos (flu) ’ 
sinh(Art) 
cos (fiu) ’ 


uj 5 (u) = oj 6 (u ) = 1, 
to 7 (u) = uj$(u) = tan (/zw), 


( 2 . 20 ) 


where Aare parameters. The RE in component forms is equivalent to the following twelve equations: 


f u 2 v 7 (x 2 y 2 + x 3 y 3 ) + u 5 Vi(x 3 y 2 - x 2 y 3 ) = 0 , 

\ u 7 v 2 {x 2 y 2 + x 3 y 3 ) + UiV 5 (x 3 y 2 - x 2 y 3 ) = 0 , 

f u 7 v 1 (x 1 y 1 - x 4 y 4 ) + uru 7 (x 4 yi - xiy 4 ) = 0, 

\ u 2 v b (xiyi - x 4 y4) + u 5 v 2 (x4yi - xiy 4 ) = 0, 

! u 5 v 2 xiy 2 - uiv 7 xiy 3 + u 2 v 5 x 4 y 2 
u 3 v 2 xiy 3 - u\v 7 xiy 2 + u 2 v 5 X4y 3 
u 5 v 2 X4y2 - uiv 7 X4y 3 + u 2 v 5 xiy 2 
u 5 v 2 X4y 3 - uiv 7 X4y2 + u 2 v 5 xiy 3 


- uiv 2 )x 2 yi + (u 5 v 7 - u 7 v 5 )x 3 yi = 0 , 

- uiv 2 )x 3 yi + (u 3 v 7 - u 7 v 5 )x 2 y 1 = 0, , 

- uiv 2 )x 2 y4 + (u 5 v 7 - u 7 v 5 )x 3 y 4 = 0 , 

- uiv 2 )x 3 y4 + {u 5 v 7 - u 7 v 5 )x 2 y4 = 0 , 


- u 7 vi X4y 3 + (u 2 vi 

- U 7 V\X4y2 + (u 2 Vi 

- u 7 viXiy 3 + (u 2 v 1 

- u 7 viXiy 2 + {u 2 vi 
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! {uiVi + u 2 v 2 )x 4 y 2 + {u 7 v 7 - u 5 v 5 )x 4 y 2 
(uiV! + u 2 v 2 )x 1 y 3 + ( u 7 v 7 - U 5 V 5 )x 4 y 3 
(uiVx + u 2 v 2 )x 4 y 2 + ( u 7 v 7 - u 5 v 5 )x!y 2 
(uiV! + u 2 v 2 )x 4 y 3 + ( u 7 v 7 - u 5 v 5 )xiy 3 


u 5 vix 2 yx + uiv 5 x 2 y 4 + u 2 v 7 x 3 yi - u 7 v 2 x 3 y 4 = 0 , 
u 3 v 4 x 3 yi + uiv 5 x 3 y 4 + u 2 v 7 x 2 yi - u 7 v 2 x 2 y 4 = 0 , , > 

u 3 vix 2 y 4 + uiv 3 x 2 yi + u 2 v 7 x 3 y 4 - u 7 v 2 x 3 yi = 0 , 
u 3 v\x 3 y 4 + u 3 v 3 x 3 yi + u 2 v 7 x 2 y 4 - u 7 v 2 x 2 yi = 0 . 


We find that under R -transformation of interchanging to 2 {u) and u> 3 (u) in ( 2 . 2 C), one can perform a 
//-transformation as follows 


a 2 (u) = - a 2 (u ), or a 3 {u) = -a 3 {u), 
to keep the system of equations (C) invariant. 


( 2 . 21 ) 


The existence of nontrivial solution implies that there exists relation 

u 7 v 2 _ u 4 v 3 
u 2 v 7 U 5 Vi ’ 

which requires A = dziy. Thus we should consider two different //-matrices, 


= u 4 (u) = 1, 

lu 2 (u) = — u> 3 {u ) = i tanu, 
u 5 (u) = u 6 (u) = 1 , 
u> 7 (u) = u>s(u) = tan u, 


( 2 . 22 ) 


and 


f Wi(lt) = u 4 (u) = 1, 

W 2 {u) = —ui 3 {u) = —i tanu, 
uj 5 (u) =u> 6 (u) = 1 , 
u ) 7 (u) = u>s (u) = tan u. 


( 2 . 23 ) 


They are in fact related each other by an exchange lu 2 (u) <-> u 3 ( u). L et us give solution iv_(u) directly 
because the calculation procedure has nothing new. For R(u) in ([ 2 . 22 |), we have 


while for R(u ) in 


K (u) — ( Mi(l + v i sin2rt) iy 2(1 + v 2 cos2 u) sin2rt 
y A*2(1 — ^2 cos2u) sin2u y\(l — v 4 sin2u) 

, using //-transformation (}2.2l|), we have 


( 2 . 24 ) 


K_ ( u) 


/ri(l + v\ sin2u) 

— y 2 (l — v 2 cos2w) sin2w 


iy 2(1 + v 2 cos2w) sin2w 
yi(l — v\ sin2w) 


( 2 . 25 ) 


where yi, y 2 , v \, v 2 are free parameters. 


AV-Matrix for Six-Vertex Model 


The general six-vertex //-matrix takes the form 


R(u) = 


( wi(w) 

0 

0 

V 0 


0 

u 2 (u) 
UlQ (u) 

0 


0 

w 5 (m) 

u 3 (u) 

0 


0 
0 
0 

u> 4 (u) J 


( 3 . 1 ) 


By setting u 7 $ = 0 and v 7t g = 0 in Eqs.(A), we write down the reflection equations for 6V type 
//-matrix in component forms: 


(u 4 - ui)v 2 x 2 y 2 = 0 , 
(«4 - ux)v 3 x 3 y 3 = 0 , 


(D.l) 









uiv 5 (x 2 y 3 - 2:32/2) = 0, 
u 4 v 6 (x 2 y 3 - x 3 y 2 ) = 0, 


(D. 2 ) 


f {u 2 V3 - u 3 v 2 )x 4 y 4 + u 5 vi(x 3 y 2 - x 2 y 3 ) = 0, 

\ (2x2X3 - 2x3X2)2:42/1 + u 6 v 4 (x 3 y 2 - x 2 y 3 ) = 0 , 

f u 3 v 6 x 1 y 1 - u 3 v b x 4 y 4 - u 6 v 3 Xiy 4 + u 3 v 3 x 4 yi + u 3 {v 4 - Vi)x 3 y 2 = 0 , 

\ u^Xiyx - u 2 v 5 x 4 y 4 - u G v 2 x 1 y 4 + 1x5X22:42/1 + u 2 {v 4 - Vi)x 2 y 3 = 0 , 


{ 1x6X22:12/2 + U 2 v b x 4 y 2 + (u 2 v 1 - Uiv 2 )x 2 y 1 = 0, 

1x6X32:12/3 + 2x3^5X42/3 + (2x3X1 - u 1 v 3 )x 3 yi = 0, 

2x5X2X42/2 + u 2 v 6 Xiy 2 + (ix 2 x 4 - u 4 v 2 )x 2 y 4 = 0, 

X5X3X42/3 + 2X 3 X 6 Xi 2/ 3 + (ix 3 x 4 - x 4 x 3 )x 3 2/4 = 0, 

! (lXlXi - M 3 X2 )xi 2/2 - 2X5X5X42/2 - 2X5X1X22/1 + 2X1X52:22/4 = 0, 

(2X1X1 - 2X2X 3 )xi 2/3 - 2X5X5X42/3 - 2X5X1X32/1 + 211X52:32/4 = 0, 

(1X4X4 - U 3 v 2 )x 4 y 2 - u 6 v e xiy 2 - u & v 4 x 2 y 4 + 2x4X6X22/1 = 0, 

(1x4x4 - 1x2X3)3:42/3 - X6X6X12/3 - X6X4X32/4 + X4X6X31/1 = 0. 


(D. 5 ) 


(D.6) 


From Ref. |^T| , we know that the 6V type solutions of YBE have the same solution-transformation as 
that for 8V type solutions except for the symmetries of weight functions and of the interchanging indices 
related to 0x7(it) and uj 3 {u). Now the two symmetries of weight functions are 

Ul 2 (ix) = SOX 2 (lt), U>3 (ix) = S" 1 0X3(21) ( 3 . 2 ) 


and 

015(21) = e cu u> 3 {u ), a>6(ix) = e~ cu coe(u) ( 3 . 3 ) 

where s,c are two nonzero constants. In fact, we find that transformation (|]^) has no effect on the 
system of equations (D), and if making ^'-transformation 

01(21) = e c “ai(ix), 04(12) = e~ cu a 4 (u) ( 3 - 4 ) 


the new A'_(«) is still a solution to RE for the new A-matrix obtained from A-transformation ( 3 . 3 ). Due 
to these symmetries, we will consider the gauge A-matrices as in 8V model. They are also classified into 
two classes- the Baxter type 

l( 21) = w 4 (ix) =- \- h ' , 


and the free-Fermion type 


“iW-^w - —55 W 
^(1X)=- 3 (2X) = S 

0X5(11) = tX 6 (ix) = 1, 


, \ sin(ix + h) 

^ = sin h \ 

, x sin(—ix + h) 

^ = sin/x » 

“*(«) = **(«) = -» 
W 5 (ix) = Uq (ix) = 1. 


( 3 . 5 ) 


( 3 . 6 ) 


For Baxter-type, the general solution to RE was given in 

A_(ix) = 


A sin(a — ix) /x sin( 2 ix) 
v sin( 2 ix) A sin(a + ix) 


( 3 . 7 ) 


which has four free parameters A, a, /x, and v. 


While for free-fermion type, since 021(21) 7^ 0x4(it), one can immediately see that o 2 (ix) = 0 and 
03(21) = 0 from (D.l). In other words, the RE for the free-Fermion type 6V models only has diagonal 
solution, 

sin(a — ix) 

0 


I\-(u) = 


9 


0 

sin(a + ix) 


( 3 . 8 ) 










In addition, if setting cos/i = 0 in (3.6), we have symmetric A matrix of free-Fermion type as follows 


R(u) = 


( cos u 0 0 0 

0 sinrt 1 0 

0 1 sinrt 0 

^ 0 0 0 cosit J 


(3.9) 


Just as discussed in Remark 2.3.2, this A-matrix shares the same AT_-matrix in (BJ) with 6V Baxter-type. 


So, up to AT-transformation (3.4), we obtain all general solutions ( p.7| ) and ( p.8[ ) to reflection equation 
in six-vertex case. 


4 A"_-Matrices to Seven-Vertex Model 


If setting u)g(u) = 0 in eight-vertex A-matrix (E3>* we get seven-vertex one 


R(u) 


( u>i(u) 0 

0 u) 2 (u) 

0 ojq{u) 

\ 0 0 


0 u 7 (u) \ 

lu 5 (u) 0 
oj 3 (u) 0 
0 (04(14) ) 


(4.1) 


The classification of solutions to the coloured 7V-type YBE is given recently in ^6| . Due to less symme¬ 
tries, the A-matrices show a much more different properties from that of both eight-vertex and six-vertex 
models. At this point, we expect that the corresponding reflection equation reveals new features as well. 


First of all, let us study symmetries of reflection equation as do previously for other cases. After 
removing the terms containing u 3 and v 3 in the system of equations (A), we find that there still exists 
A'-transformat ion ( |2.4| ) under A-transformation ( |2.3| )(note that ujs is absent!). 

In ^(|, an additional relation oj 3 (u) / u> 3 (u) = e cu is given, where c is a constant. When c ^ 0, there 
have only trivial A'_-matrices. The case of c = 0 or 0 ) 5 ( 11 ) = u>e(u) is further classified into three different 
types: Baxter type, free-fermion type I, II, which will be discussed in the following sections. 


4.1 Baxter Type 


The parametrization of the A-matrix is as follows 


/ \ / \ sin(it + h) 

Wl (u) = W 4 («) = - 

**(„) = «,*(„) = -gs, 

U) 5 (u) = lo 6 (u) = 1 , 
oj 7 (u) = sin(it + h) shut. 


(4.2) 


Substituting (4.2) into (A), we solve these equations case by case. 


Case 4-1-1: Diagonal Solution. It can be seen from (A.l) that if 02 ( 14 ) = 0 then 03 ( 14 ) = 0 (Note that 
us = 0 = us in (A)). We obtain the diagonal solution: 


A'_(u) = 


sin (a — 14 ) 


sin(a + 14 ) 


(4.3) 


Case 4-1-2: Skew-diagonal solution. Let a\(u) = 0, it requires 04 ( 14 ) = 0 from (A.4). We only need 
to consider two equations: 

U 7 i’ 3 x 3 y 3 + u 1 v 5 (x 2 y 3 - x 3 y 2 ) = 0 , 
u 3 v 7 x 3 y 3 + u 3 v 1 (x 2 y 3 - x 3 y 2 ) = 0 . 


(4.4) 
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Solving Eqs.(4.4), we have two A'_-matrices, 


(! t') 


where p(u) = (A + cos2w)/2 and A is a free parameter. 

Case 1.3: 03 (u) = 0. One can get from case 4.1.1 that 

ai(tt) = p(u) sin(a — it), a 4 (u) = p(u) sin(a + u ) 


and find 

so the AC-matrix is 

where p, v are parameters. 


o 2 {u)/p{u) = p sin(2u), 


K_{u) = 


v sin(a — u) p sin( 2 u) 

0 v sin(o; + u) 


( 4 . 5 ) 


Case 4 - 1 - 4 - Combining the results obtained above, one can easily write down the general A'-matrix 
as follows 

.. ’ (4-6) 


K = ( v sin (a - u) pp(u) sin(2u) A 
~~ ' ' p sin(2u) v sin(a + u) j 


In summary, we can regard ( |4.5| ) and (4.6) as the most general reflection matrices, because others 
can be obtained by assigning special values to free parameters. Furthermore, comparing (4.3) and (4.5), 
we see that in the case of 7V type, 02 (u) = 0 implies 03 (u) = 0, but the reverse does not hold, this is 
different from the case of 8 V type. 


4.2 Free-Fermion Type I 


This kind of A-matrix reads 


, s sin(u + h) 
^ = sin h ’ 

, sin (—u + h) 
^ = sin h ’ 

u 5 (u) = u/ 6 (u) = 1, 


( 4 . 7 ) 


If cosh = 0, or wi(it) = u) 4 (u), we can make the similar calculation as do for the Baxter-type, for both 
have the same symmetries. The result is almost the same as that given in ( 0 ) and ( f4.5| ) but with 
p(u) = (A + cos 2 m). 

When u>\{u) / (it), it force 03 ( 1 /) = 0 from the second equation of (A.3). The RE reduces to the 

following equations 


r u 2 v 5 (xxyi - X4J/4) + u 5 v 2 (x 4 yx - *12/4) = 0, 

\ u 7 vxxiy 1 - u 7 v xx 4 y 4 + u 4 v 7 x 4 yx - u 4 v 7 xxy 4 + (u 4 - u 4 )v 2 x 2 y 2 = 0, 

! U5V 2 xiy2 + u 2 v 3 x 4 y 2 + {u 2 vi - u\v 2 )x 2 yx = 0 , 

U5V 2 X4y 2 + U 2 V 5 Xiy 2 + ( m 2^4 - U 4 V 2 )x 2 y 4 = 0 

(itiVi - ■U 3 v 2 )x 1 y 2 - u 5 v 5 x 4 y 2 - u 5 vxx 2 yx + uxv 5 x 2 y 4 = 0 , 

{u 4 v 4 - u 3 v 2 )x 4 y 2 - u 5 v 5 xiy 2 - u 5 v 4 x 2 y 4 + u 4 v 5 x 2 y 4 = 0 . 


Case 4-2.1: a 2 {u) = 0. The two equations of (E.l) are not compatible with each other due to the 
symmetry between u>x(u) and u> 4 (u) being broken. Therefore there exists no diagonal solution for this 
type of A-matrix. 
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Case 4-2.2: a\(u) = 0. One can deduce from (E.l) that 444 ( 14 ) = 0 and a 2 {u) = 0 . This is a trivial 
case. 

Case 4-2.3: General solution. From the second equation of (E.l) and (E.2), one can get the same 
result as that in Baxter type 

ai(rt) = sin(a — u), 444 ( 14 ) = sin(a + u), a 2 (u) = /i sin2«. 

Substituting them into the second equation of (E.l), one find fi = ±1. So, we have only one solution 


sin(cx — u) ± sin(2w) 
0 sin(a + u) 


K-(u) = 

which also shows that 03 (44) = 0 does not imply a 2 (it) = 0 . 


(4.8) 


4.3 Free-Fermion Type II 

In this case, the elements of the P-matrix take the following forms 

{ u>i[u) = 074(44) = cosh it, 

U 2 {u) = -ui 3 (u) = sinhu, 

u 5 (u) = u 6 (u) = 1 , 

W7 (u) = u. 

For the sake of brevity, we simply give the result. Note that the nontriviality requires 444 ( 14 ) = ± 04 ( 14 ). 
If ai(u) = 04 ( 14 ), we get 

*-(«)=(“ t ‘ ai ‘ hU ). (4.10) 

while if 01 ( 14 ) = — 04 ( 14 ), we have 

*-(»>=(„ ^ u ). (nil 

In addition, according to the discussion in section 2.4, the solutions keep invariant under exchange of 
u> 2 <-> u >3 since 03 ( 14 ) = 0. 


5 Construction of Boundary Hamiltonian 


In this section, we will discuss the Hamiltonians for the systems described by the P-matrices and K- 
matrices obtained in the previous sections. The 6V (Baxter type and free-fermion type) and 8V (Baxter 
type and free-fermion type I) are included in Sklyanin’s formalism. While for 7V models, both Baxter 
type and free-fermion type I P-matrices has only regularity, P-symmetry, unitarity and crossing-unitarity 
symmetries, their K + (i4)-matrices are obtained by (1.5). However, all of these cases has the same defini¬ 
tion of transfer matrix |q, Il2l, and we can construct their Hamiltonians in a unified way. 


If A±(0) oc id , trK + ( 0) ^ 0, the Hamiltonian for the open systems is defined asj 


N~ 1 1 1 

n = Y. + 2 A - 1 (°) K '~ (°) + 

3 =1 

where two-site Hamiltonian TLjj+i is given by 


tr 0 I< + (0 )H N q 
trK+( 0) 


(5.1) 


^.7 J + 1 1 (^) lit —0 (ju A. 7 ,i +1 C 1 ^) I u —0 Pj,j+ 1 - ( 5 - 2 ) 

All the Baxter type models in two-component systems belong to this case. The boundary Hamiltonian 
of 6V and 8V Baxter type can be found in |l3j, [li| [l^] . 
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Case 5.1: Baxter type 7V with crossing parameter 77 = h. From K_(u) in (4.6) and relation (111), we 
find that 

K + (u) = K_(-u- h-,-a + ,n+,i/ +) A+). (5.3) 

According to Eqn.(|5.l|), the Hamiltonian is 


JV-l 


Ft — ——7—— y 7~tj 7+1 — A—u\ + B—C rl + C—C 
4 sin h ^ 


i=i 


—A + <j~ n + B + a + C+cr^Y, 


(5.4) 


where 


ftjj+i = (2 + sin 2 h)aj crj +1 + (2 - sin 2 /i)crjcrj +1 + % sin 2 /i((j|(Tj +1 + crjcr^ +1 ) + 2 cos/iu|cr| +1 , (5.5) 

A± = ^ cota±, B ± = (1 + A±)M± , C'+ - M . (5.6) 

z z sina± z/± sina± 


However, if trA"+(0) = 0, just as pointed out in Refs.JI^, |l 6 ), there will have no well-defined Hamil¬ 
tonian from the first derivative of the transfer matrix as done in ( |5. l[ ) . But if 

tr 0 K+ (0)7"tjv O = A -id, (5.7) 

where A is a constant, we can still derive the well-defined local Hamiltonian from the second derivative 
of transfer matrix as follows 

P'fnl JV_1 1 1 

« S 4icTi4) = g H " + ' + 2 A ' r,(0) ^- (0) (5 - 8) 

+ 2(C + 2A) {t r o{K+ {0 )Gno) +2tr 0 (K' + (0 )Hno) +tr 0 (K+ (0 )FL 2 NQ )}, 

where 


C = trK' + ( 0), (5.9) 

n — p ^ ^3,3+l( u ) I (r in, 

*3,3+1 = fj,j +1 -- 1 «= 0 - (O.IU) 

The following discussions show that all the boundary conditions corresponding to the free-Fermion type 
A-matrix belong to this case. We argue that it is a common property for all free-Fermion models. 


Case 5.2: Free-Fermion type-I 8V with crossing parameter 77 = I. Here I is the complete elliptic 
integral of the first kind of modulus k. For general boundary condition described by A'_(it) in (2.17), we 
have 


K+(u) = Kti-u-I) 


k ,2 F+(u) snit + E+(u) cnit dnit 2 p+k' 2 snit cnit dnit 

2 ep+k 12 snit cnit dnit k l2 F+(u) sn u — E+(u) cnit dnit 


where 


A+(it) = cf dn it + 


+ j„ 2 „. , k({l-ekG)c+ + HcJ) 


eG — k 


cn u , 


„ ^ , 2 , m ~ ekG)c+ - k /2 Hc+) 2 

A+ (it) = c 2 dn it -\ -—- 7 -cn it. 


eG-k 


From Eq.(5.S), the Hamiltonian reads 


AT —1 


Ft — 'y ' Ftjj+i + A_erJ + H_(cr^ + ecr 1 ) + A+cr]^- + H+(cr)^ + eu N ), 

3=1 


(5.11) 


(5.12) 

(5.13) 


(5.14) 


13 













where 


(5.15) 


JCj- 

Kjj+i = y (°7 + a j+ 1 ) + 


G + k 


G-k 


2 1+ 2 


and 


^ 4 - = c 2 / c i , B- = 2 fi /c x , 

A+ = k' 2 (HF + (0) - A+(0))/2(fc , 2 A + (0) F HE+( 0)), 
A+ = fc , 2 (G' + e/fc)/r+/(A; , 2 F + (0) + HE + ( 0)). 


For the diagonal A'-matrix (2.12), we have 


JV -1 


,_ aV 

B = ^2 ^ 1 , 1 + 1 + ~( a l ~ a N)i 
1=1 


(5.16) 


which has been discussed in pi and is a special case of (5.14). 


Case 5.3: Symmetric free-Fermion type-I 8V. If considering A-matrix ( [2.1H| ) and the corresponding 
A'-matrix(2.9), we get the following Hamiltonian 


N-l 


1 + k 


H = 


1 - k 


1 u !+1 


- V rr V 

J j <J j+ 1 


1=1 


where 


A—(j x F B—(j// + CActj H- -f- C+cTjy 


. cn a± dna± , cna + , 9 . 

A_ = —-, A+ = -- f - (1 — k 2 sna+), 

2 sna-t 2 sna + dna+ 

= /i±(A ± + 1) ^ _ ^±(A± - 1) 


(5.17) 


sna;± 


sna± 


Case 5.4: Free-Fermion type 6V with ij = 7 t/ 2 . For A-matrix in dU) and the general A'_-matrix 
(3.8), if setting 

K + (u) = K*_{-u — 7 t/ 2 ; tt/ 2 — a+ — h) (5.18) 


we have 


N -1 

n = —^ E ( a i CT "+i + fyj+i + cos/l H + a l+i)) 

3 = 1 

— cota_ai — cota+af^. 


(5.19) 


Case 5.5: Symmetric free-Fermion type 6V. If considering R(u) in (3.9) together with the general 
K --matrix (3.7), we can set 

K + {u) = K l _(-u - tt/2; -a + ,p + ,u + ), 


thus the Hamiltonian is 


N-l 


n = y i - « i * 

1=1 

A—cr x F B-<j x F G-g x A+fj^ F B+cr~p] -f- C+Cpj 


where 


A± = cota±, B± = 


2 /^± 

sina± 


C± = 


2v± 

sina± 


(5.20) 


(5.21) 


Case 5.6: Free-Fermion type I 7V with crossing parameter r] = -k/2. From A'_-matrix in (^(^) and 
relation (HI), we get 

K + (u) = K-(—u — 7 t/ 2 ; a+ — h + tt/2, fx+), (5.22) 
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and the Hamiltonian is 


where 


H 


2^ 51 { cosh K + CT I+i) + 2a J a J+i + *(°i °i+i + °i°j+i)} 

i=i 

—+ A_cr^ — + A+crjy, 


/i_,/u+ = ±1, A± 


cotci'±/2, B± 


sina± 


(5.23) 


Case 5.7: Symmetric free-Fermion type I TV. From AY-matrix in (4.6) with p(u) = A + cos2 u and 
relation (|1.6|), we have 

K+(u) = K_(-u - n/2-, -a+, n+,v+, \ + ), (5.24) 


and 


where 


n 


N -1 


E^ x 


3 U 3 +1 


3 =1 


+ 2^1 +1 + ^i)} 


A_u* T B-di + (7_(J 1 + A+crjy -f- Clf-u^y 


A± 


cota± (Id - A±)/U± _ /^± 

^ ? B± : ? C'± : • 

2 v± sina-t v± sina± 


(5.25) 


It should be pointed out that for the free-Fermion type II of both 7V and 8V models which have no 
crossing-unitarity symmetry, how to prove their integrability and to obtain the corresponding Hamilto¬ 
nians in the case of open boundary condition is an open problem. 


6 Remarks and Discussions 


In this paper we find that symmetries play an important role in solving the reflection equation. For 
any non-standard A-matrix which is obtained by applying /^-transformation to the standard one, then 
the corresponding reflection matrix can be obtained by making A'-transformation to that for standard 
A-matrix. 

Moreover, all solutions given above indicate that the number of free parameters appeared in AT-matrix 
is determined by symmetries of A-matrix. The A-matrices with different forms but the same symmetries 
share the same AT_(u) matrix. The free-Fermion type A-matrix with u>\{u) = u> 4 (u) is just in this case. 
It has the same form A'_(w) as in Baxter type. Also we note that, different from that for six- and eight- 
vertex cases, the elements a 2 {u),a^[u) of I\-(u ) in seven-vertex case have no interchanging-symmetry 
resulting from the symmetry between lot(u) and of A-matrix being broken. 

It is also interesting to note that while constructing Hamiltonian, all reflection matrices for free- 
Fermion A-matrices have property of trK+(0) = 0. We argue that it is a typical property for all 
free-Fermion models. So the local Hamiltonian for such system are obtained from second derivative of 
the transfer matrix. 

We are sure that our procedure to find solutions of the reflection equation can be applied to high-spin 
models, though the calculation may be much more involved in this case. With the solutions given in this 
paper, we can use the Bethe ansatz method to study the physical properties of open spin chains. 

Furthermore, recently much attention has been directed to the Yang-Baxter equation with dynamical 
parameters |27], |28[. How to construct the corresponding reflection equation and to seek its solution is 
an open problem. We wish to discuss some related problems using the method and procedure given in 
this paper. 
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